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ABSTRACT 

We present Hubble Space Telescope ultraviolet spectroscopy of the white dwarfs 
PG 0843+516, PG 1015+161, SDSS 1228+1040, and GALEX 1931+0117, which ac- 
crete circumstellar planetary debris formed from the destruction of asteroids. Com- 
bined with optical data, a minimum of five and a maximum of eleven different metals 
are detected in their photospheres. With metal sinking time scales of only a few days, 
these stars are in accretion/diffusion equilibrium, and the photospheric abundances 
closely reflect those of the circumstellar material. We find C/Si ratios that are consis- 
tent with that of the bulk Earth, corroborating the rocky nature of the debris. Their 
C/0 values are also very similar to those of bulk Earth, implying that the planetary 
debris is dominated by Mg and Fe silicates. The abundances found for the debris 
at the four white dwarfs show substantial diversity, comparable at least to that seen 
across different meteorite classes in the solar system. PG 0843+516 exhibits significant 
over-abundances of Fe and Ni, as well as of S and Cr, which suggests the accretion 
of material that has undergone melting, and possibly differentiation. PG 1015+161 
stands out by having the lowest Si abundance relative to all other detected elements. 
The Al/Ca ratio determined for the planetary debris around different white dwarfs is 
remarkably similar. This is analogous to the nearly constant abundance ratio of these 
two refractory lithophile elements found among most bodies in the solar system. 

Based on the detection of all major elements of the circumstellar debris, we cal- 
culate accretion rates of ~ 1.7 x lO^gs"^ to ~ 1.5 x lO^gs"^. Finally, we detect 
additional circumstellar absorption in the Si iv 1394,1403 A doublet in PG 0843+516 
and SDSS 1228+1040, reminiscent to similar high-ionisation lines seen in the HST 
spectra of white dwarfs in cataclysmic variables. We suspect that these lines originate 
in hot gas close to the white dwarf, well within the sublimation radius. 

Key words: Stars: individual: PG 0843+516, PG 1015+161, 

SDSSJ122859.93+104032.9, GALEX J193156.8+011745 - white dwarfs - cir- 
cumstellar matter - planetary systems 



1 INTRODUCTION 

Most of our current insight into the interior structure of exo- 
planets is derived from the bulk density of transiting planets 
(e.g. Valencia et al. 2 010 ). and transit spectroscopy provides 
some informati on on the chemical co mposition of their at- 
mospheres (e.g. iGrillmair et al.ll200^ '). More detailed inves- 
tigations of the chemistry of exo-planetary systems around 
main-sequence host stars are beyond the reach of present 
observ ational instrumentation. However, IZuckerman et all 
(|2007l ) demonstrated in a pioneering paper that the pho- 
tospheric abundances of polluted white dwarfs can be used 
to infer the bulk abundances of the planetary debris mate- 



rial detected around the white dwarf GD 362, and showed 
that the composition of this material is broadly comparable 
to that of the Earth-Moon system. 

The strong surface gravity of white dwarfs implies that 
metals will sink out of the photosphere on time scales that 
are orders of magnitude shorter than their cooling ages, and 
therefore white dwarfs are expe cted to have either pure hy - 
drogen or helium atmospheres (|Fontaine fc Michaudl[l97a \ 
Exceptions to this rule are only hot (Tcff ^ 25 000 K) white 
dwarfs where radiative levitati on can support some heavy el- 
ements in the photosphere (e.g. lChaver et aDl995l ). and cool 
{Tcs < 10 000 K) white dwarfs where convection may dredg e 
up core material (jKoester et al.|[T982l : iFontaine et al.lll984 ). 
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Yet white dwarfs with metal-contami nated atmosphere s 
have been known for nearly a century (|van MaanerJ 191^ . 
and a ccretion from th e interstellar medium (e.g. iKoested 
1 19761 : IWesemaell 1 19791 : iDupuis et al] Il993l ) has been the 
most widely accepted scenario, despite a number of funda- 
ment a l problems (e.g.lAannestad et al.lll993l : iFriedrich et alj 
l2004l : iFarihi et al. I l2010ah . However, the rapidly growing 
number of white d warfs that ar e accr e ting from c i rcum- 
stellar discs ( e.g. iBecklin et alJ 120 05: Kilic et alJ l2005l: 



Gansicke et al 



200e 



20081 : IVennes et al 



von Hippel et al.. 2007: Farihi et al.l 



2010l : iDufour et al.1 \2Qm unambigu- 
ously demonstrates that debris from the tidal disrup- 
tion of main-belt analogue ast eroids or minor planets 
l|Graham et"al] Il990l: I Jural 120031 '). or Kuiper-belt like ob- 



iects (iBonsor et al.ll201l|) . likely perturbed by unse en plan- 



ets (|Debes fc SigurdssonI I2OO2I : iDebes et all \20li ). is the 



most likely origin of photospheric metals in many, if not 
most polluted white dwarfs. 

Because of the need for high-resolution, high-quality 
spectroscopy, detailed abundance s tudies have so far been 
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dance and white dwarf temperature, metal lines are stronger 
in a helium-dominated (DB) atmosphere than in a hydrogen- 
dominated (DA) atmosphere, as the opacity of helium is 
much lower than that of hydrogen. Therefore, the small sam- 
ple of well-studied metal polluted white dwarfs is heavily 
biased towards DB white dwarfs, which have diffusion time 
scales of ~ 10^ — 10^ yr. These long diffusion time scales 
introduce a significant caveat in the interpretation, as the 
abundances of the circumstellar debris may substantially 
differ from those in the white dwarf ph otosphere if th e accre- 
tion rate varies on shorter time scales (|Koesterll2009l ). While 
the life times of the debris discs arc subject to larg e uncer- 
taint ies, there are theoret ical (j Rafikov 2011; Mctz ger et al.l 
I2OI2I ) and observational l|Girven et al.l 120121 . Farihi et al. 
2012 in press) arguments that suggest that the accretion 
rates onto the white dwarfs may vary significantly over pe- 
riods that are short compared to the diffusion time scales. 
In fact, some of th e most heavily polluted white dwarfs have 
no infrared excess (IFarihi et al ] |2009l : iKlein et ahlbOllT ) , and 
may have accret ed all the circumst e llar debris a few dif fusion 
time scales ago l|Farihi et al.ll2009l : iGirven et al.ll2012l 'l. 

We are currently carrying out an ultraviolet spectro- 
scopic survey of young DA white dwarfs that have cooling 
ages of 20 to 200 Myr, metal sinking time scales of a few 
days, and are hence guaranteed to be in accretion-diffusion 
equilibrium. The aim of this survey is to determine the frac- 
tion of white dwarfs that are presently accreting planetary 
debris, and to determine accurate abundances for a subset. 
Here we present the analysis of four heavily polluted white 
dwarfs that are known to also host planetary debris discs. 



2 OBSERVATIONS 

The targets for our ongoing far-ultraviolet spectroscopic 
survey of young and correspondingly warm (17 000 K < 
Tcfi < 25 000 K) DA white d w arfs w ere dr awn from the 
compi lations of iLiebert et al.l l|2005h and iKoester et al.l 
l|2009l ), supplemented with a few recent discoveries (e.g. 



iGansicke et al.l I2OO6I : IVennes et al.l I2OI0I ) . Our sample also 
includes a small number of post-common envelope binaries 
(PCEBs) in which the white dwarf accretes from the wind of 
the M-dwa rf companio n. Th ese s ystems were selected from 
ISchreiber fc Gansick^ (120031 ) and IFarihi et all ([2010b) with 
the same cut on white dwarf temperature and cooling age. 
Under the assumption that the M-dwarfs have a solar-like 
composition, the white dwarfs in PCEBs serve as "abun- 
dance standards" for our abundances analyses and diffusion 
calculations. 

2.1 HST/COS spectroscopy 

PG0843-f-516, PG1015-H61, and 

GALEX J193156.8-F011745 (henceforth 
GALEX 1931-1-0117) were observed as part of our snapshot 
survey, with exposure times of 1420 s, 1424 s, and 800 s, 
respectively. We used the G130M grating with a central 
wavelength of 1291 A, which covers the wavelength range 
1130 - 1435 A, with a gap at 1278 - 1288 A due to the 
space between the two detector segments. To mitigate the 
fixed pattern noise that is affecting the COS far-ultraviolet 
detector, we split the exposure time equally between 
two FP-POS positions (lfc4, the limited duration of the 
snapshot visits did not allow to use the full set of four 
different FP-POS positions). 

We also report COS observations of three PCEBs 
observed within this snapshot survey, that will be 
used as "abunda nces standards": GD448 (HRCam , 
iMaxted et al]|l998l ). CD 245 (MS Peg, 'Schmidt et al."l995l), 
and PC 2257-^162 (KUV22573-f 1613, Wachtcr ct al. 200:|), 
with exposure times of 900 s, 600 s, and 1070 s, respectively. 

SDSS J122859.93+104032.9 (henceforth 
SDSS 1228-1-1040) was observed in Cycle 17 as part of 
a regular Guest Observer programme. We obtained two 
sets of spectroscopy with the G130M grating with central 
wavelengths of 1291 A and 1327 A, and both observations 
were again split among two FP-POS positions (lfc4). In 
addition, we obtained G160M spectroscopy with central 
wavelengths of 1577 A and 1623 A. The total exposure time 
of the G130M and G160M observations were 2821s and 
4899 s, respectively, seamlessly covering the wavelength 
range 1130 - 1795 A. 

The data retrieved from the HST archive were pro- 
cessed and calibrated with CALCOS 2.15.6. The COS spec- 
tra of the four white dwarfs shown in Fig.[T] reveal the 
broad Lya profile typical of DA white dwarfs, plus a mul- 
titude of narrow absorption lines from a range of metals. 
The peak signal-to-noise ratio (S/N) of the COS spectra 
is reached in a line-free region near 1320 A, and ranges 
from ~ 25 for PG0843-H516 and PC 1015-^161 to ~ 40 for 
SDSS 1228-H040 and GALEX 1931+0117. However, these 
values only include photon count statistics, and do not ac- 
count for the residual fixed-pattern noise related to the use 
of only two FP-POS positions. The resolving power of the 
COS spectra, as measured from on-orbit data ranges from 
- 15 000 at 1150 A to ~ 20 000 at 1430 A. 



2.2 Optical observations 

The wavelength spanned by our COS observations does not 
cover any strong line of either Ca (traditionally the most 
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Figure 1. COS/G130M spectra of four white dwarfs known to iiave circumstellar discs, scaled to a peak flux of unity, offset by 1.4 units, 
and sorted from top to bottom by increasing metal abundances. For warm white dwarfs with pure hydrogen atmospheres, the broad 
Lya line is the only spectral feature in this wavelength range. These four stars accrete from the circumstellar debris, and their spectra 
are riddled with absorption lines of C, O, Al, Si, P, S, Cr, Fe, and Ni. In addition, Mg and Ca can be detected in their optical spectra 
(Fig.©. 



important tracer of metal pollution in white dwarfs, and 
an important refractory element) or Mg (one of the major 
constituents of rocky material in the solar system, includ- 
ing the Earth). Ground-based abundance studies using the 
Ca II H/K doublet and the Mg ii 4482 A hne ar e already pub- 
lished for GALE X 1931-H0117 (IVennes et al.1 12010| . l2011bl : 
iMelis et al.ll201ll '). Two short (lOmin) VLT/UVES spectra 
of PG 1015-1-161 w-ere o btained as part of the SPY p roject 
ijNapiwotzki et al.lliooil 'l. which iKoester et al.l l|2005f l anal- 
ysed to determine the Ca abundance of PG 1015-1-161 
(Sect. l4^ . Here we use the same spectra to determine in 
addition the abundance of Mg. 

We observed PG 0843+516 for a total of 2 h on the 
WHT using ISIS with the R600B grating and a 1" slit, cov- 
ering the Ca and Mg lines at a resolving power of ~ 2500 
and a S/N of s a 90. The data were reduced and calibrated 
as described in iPvrzas et af] \2Qli ). 

We also obtained a total of 9h VLT/UVES spec- 
troscopy of SDSS 1228-1-1040 between 2007 and 2009 us- 
ing the Blue390 and Blue437 setup with a 0.9" slit, cov- 
ering both the Ca and Mg features with a resolving power 
of ~ 40 000. The data were reduced in Gasgano using the 
UVES pipeline. The individual spectra were of relatively low 
S/N, and we analysed only the error- weighted average spec- 
trum, binned to 0.05 A, with S/N ~ 35. 

The optical spectra around the Call K and Mg ll 4482 A 
lines are shown in Fig. (2] We note that while most previ- 
ous studies of metal-polluted white dwarfs have focused on 
the Call H/K lines, their strength for a given abundance 
decreases strongly with increasing temperature, as Call is 
ionised to Cam. For temperatures T^s ^ 20 000 - 25 000 K, 



Table 1. Atmospheric parameters from spectroscopy 



Object 



Teff [K] logg [cgs units] 



PG 0843-1-516 = WD 0843-|-516 
optical, ^iebert_e^^l^ ^2005) 



HST, this paper 

PG 1015-1-161 = WD 1015-1-161 
optical, _Liebert_ ct al. (2005.) 
optical, J<oester_^t_aL (2009) 



23 870 ± 392 
23 095 ± 230 



HST, this paper 

SDSS J122859. 93-1-104032. 9 = 
optical, ^isenstei^_et_aL (2006) 
optical, ^ansicke_et_aL (2007) 

optical, our fit to SDSS spectrum 
HST, this paper 
adopted, this paper (Sect. 13.111 
GALEX J193156. 8-1-011745 = 
optical, ^ennes_^t_aL (2010) 
optical, Melis et al. (2011 ) 
HST, this paper 



19 540 ± 305 

19 948 ± 33 

19 200 ± 180 
WD1226-I-110 

22 125 ±136 
22 292 ±296 

22 410 ± 175 

20 565 ± 82 
20900 ± 900 

WD 1929+012 

20 890 ± 120 

23 470 ± 300 

21 200 ± 50 



7.90 ± 0.05 
8.17 ±0.06 

8.04 ±0.05 
7.925 ± 0.006 
8.22 ±0.06 

8.22 ±0.02 
8.29 ±0.05 
8.12 ±0.02 
8.19 ±0.03 
8.15 ±0.04 

7.90 ±0.03 
7.99 ±0.05 

7.91 ±0.02 



Mg II 4482 A becomes a more sensitive probe of metal pol- 
lution (e.g. iGansicke et al..,2007. : .Farihi et al.»2012ii) . 



3 ATMOSPHERE MODELS 

3.1 Effective temperature and surface gravity 

All observed HST /COS and optical spectra were analysed 
with theoret ical model atm ospheres using input physics as 
described in iKoest^ (120101') . and including the Lyma n and 
B aimer line profiles of Tremblav fc BergeronI l|2009[ ). We 
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Figure 3. The normalised COS spectra (black) of PG 0843+516 (left) and PG 1015+161 (right), along with our best-fit models (red). 
Interstellar absorption features are indicated by vertical gray dashed lines. The interstellar lines in PG 1015+161 are blue-shifted with 
respect to the photospheric features by 57kms~^, in PG 0843+516 this shift is ^ 7kms~^. Airglow of Ol can cause some contamination 
of the 1302 - 1306 A region. An illustrative airglow emission spectrum (arbitrarily scaled in flux) is shown. The strong Si IV 1394,1403 A 
doublet seen in the COS spectrum of PG 0843+516 is not of photospheric, but circumstellar origin (Sect. O. 
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Figure 4. Same as Fig.|3] but for SDSS 1228+1040 (left) and GALEX 1931+0117 (right). The interstellar linos in their spectra are are 
blue-shifted with respect to the photospheric features by 36 km s~^ and 61kms~^, respectively. The strong absorption band seen near 
1410 A in the COS spectrum of GALEX 1931+0117 is thought to be related to an autoionisation line of Sill or to a resonance feature in 
the photoionisation cross section (Sect. [372. 211 . The same features is seen, though much weaker, in SDSS 1228+1040. The Si IV 1394,1403 A 
doublet in SDSS 1228+1040 shows additional absorption, blue-shifted with respect to the photospheric features, which is of circumstellar 
origin (Sect. 
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Figure 2. The normalised optical spectra (gray) of PG 0843+516 
(WHT), PG 1015+161 (VLT/UVES, binned to 0.2 A), and 
SDSS 1228+1040 (VLT/UVES, binned to 0.05 A), along with 
the best-fit models (black). The width of the Mgll line in 
PG 0843+516 is due to the low resolution of the WHT data. In 
the spectrum of SDSS 1228+1040, the photospheric absorption of 
CallK is embedded in a double-peaked emission line from the 
gaseous debris disc, which is, however, so broad that it does not 
affect the measurement of the Ca abundance. 



used a fine grid of models spanning the range of temper- 
atures and surface gravities found for the four targets by 
previous studies (Table [1} and determined the best-fit pa- 
rameter by minimising x^, using the very good relative flux 
calibration as an additional constraint. The errors reported 
in Sect. |4] are statistical only and do not include systematic 
effects of observation, reduction, or models. More realistic 
errors can be estimated from a comparison with the other 
measurements in the literature, which used similar models, 
but optical spectra. Table [1] suggests a systematic trend for 
somewhat lower temperatures derived from the ultraviolet 
data when compared to the values based on optical spec- 
troscopy. A simil ar trend is seen for DA w hite dwarfs with 
Teft ~ 20 000 K in lLaioie fc BergeronI ()2007l '). who compared 
the efi^ective temperatures derived from optical and {Inter- 
national Ultraviolet Explorer ) ultraviolet spectroscopy. We 
carried out a range of test calculations to explore the effect 
of these systematic uncertainties in Tag and log g on the 
derived metal abundances (Sect. [331 . The abundances and 
mass fluxes do not change by more than ~ 0.1 dex, which is 
less than the typical uncertainty of our fits, and the abun- 
dance ratios vary by much less. Hence, the discussion in 
Sect. [5] and [6] is not affected by the systematic uncertainties 
in Toff and \ogg. 

Finally, to assess the possible effect that the presence of 
metals has on the effective temperature and surface gravity, 
we computed a small grid of models for the two most metal- 
polluted stars (PG 0843+516, GALEX 1931+0117), includ- 
ing metals at the abundances determined in Sect. l3.2l and 
re-fitted the HST /COS spectra. For both stars, the best-fit 
Teft and logg did not change significantly, and we there- 
fore adopted the atmospheric parameters from the pure- 
hydrogen fits for all four targets. 

3.2 Metal abundances 

The COS spectra of the four white dwarfs contain a 
multitude of absorption lines from a range of elements. 



Table 2. List of major line features used for the abundance deter- 
minations and upper limits. Because of the different wavelength 
ranges of the available spectra not all lines could be used for all 
four stars. 



Ion 



Vacuum wavelengths [A] 



Gn 
Gm 
Ni 
Oi 

Mgii 

Alii 

Aim 
Si II 



Sim 



Si IV 
Pii 



Piii 

Sii 

Siii 

Gall 

Sen 

Tim 

Viii 

Griii 

Mnll 

Mnlll 

Feii/iii 

Niii 



1334.530 
1174.930 
1199.550 
1152.150 
1239.925 
4482.383 
1670.787 
1763.869 
1379.670 
1190.416 
1250.091 
1304.370 
1348.543 
1526.707 
4132.059 
1140.546 
1155.959 
1206.555 
1312.591 
1393.775 
1149.958 
1249.830 
1543.133 
1334.813 
1250.584 
1194.041 
1169.029 
1418.773 
1298.633 
1148.465 
1136.669 
1263.611 
1162.015 
1201.118 
1174.809 



1335.660,1335.708 

1175.260,1175.590,1175.710,1175.987,1176.370 

1200.220,1200.710 

1302.170,1304.860,1306.030 

1240.395,1367.257,1367.708,1369.423, 

4482.407,4482.583 

1719.442,1724.922,1724.982,1760.106,1761.977, 

1763.952,1765.816 

1384.132,1605.766,1611.873 

1193.292,1194.500,1197.394,1246.740,1248.426, 

1250.436,1251.164,1260.422,1264.738,1265.002, 

1305.592,1309.276,1309.453,1311.256,1346.884, 

1350.072,1350.516,1350.656,1352.635,1353.721, 

1533.431,3854.758,3857.112,3863.690,4129.219, 

5042.430,5057.394,6348.864,6373.132 

1141.579,1142.285,1144.309,1144.959,1154.998, 

1156.782,1158.101,1160.252,1161.579,1206.500, 

1294.545,1296.726,1298.892,1301.149,1303.323, 

1341.458,1342.389,1365.253,1417.237 

1402.770 

1152.818,1153.995,1155.014,1156.970,1159.086, 

1452.900,1532.533,1535.923,1536.416,1542.304, 

1543.631 

1344.326 

1253.811,1259.519 
1194.433 

1169.198,1341.890,3737.965,3934.777 
1418.793 

1298.697,1298.996,1327.603 
1149.945,1149.945 

1146.342,1247.846,1252.616,1259.018,1261.865, 



1188.505,1192.316,1192.330,1197.184,1199.391, 
1233.956,1254.410 

1177.478,1179.851,1183.308,1183.863,1183.880 
many weak lines, individually recognisable 1140-1152 
1317.217,1335.201,1370.123,1381.286.1393.324,1411.065 



GALEX 1931+0117 has the richest absorption spectrum, in 
which we securely identified transitions of nine elements (C, 
O, Al, Si, P, S, Cr, Fe, Ni), and we included those metals in 
the abundance analysis of all four targets. We also include in 
the analysis N, Na, Ti, V, Mn, which have moderately strong 
transitions in the wavelength range covered by the COS ob- 
servations, but that were not detected. All metals were fully 
included in the calculation of the equation of state. 

Synthetic spectra were calculated adopting the atmo- 
spheric parameters determined in Sect . 13.11 and includ- 
ing approximately 2500 metal lines. The basic source of 
atomic line data (wavelengths, excitation energies, transi- 
tion probabilities loggf. Stark broadening constant r4) was 
V ALD (Vienna At o mic L i ne Database), which is des cribed 
in 



Piskunov et al.l (Il99^ . iRvabchikova et al.l ([l993), and 
iKupka et all (|l999l . l2000l ). The ion Sill has a large number 
of lines in the ultraviolet, and we noted a significant scatter 
in the abundances derived from different lines. Replacing the 
log gf values from VALD values with those from the NIST 
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(National Institute of Standards) database, which difTer for 
some hnes by up to 0.3 dex, leads to more consistent results. 
Nevertheless, the situation for this ion is not satisfactory 
fSect. [TT2|) . and we hav e consulted a number of original 
sources in the litera ture l|Lanz fc Artrulll985l : lNahar..l998 : 
iBautista et al.ll2009l ') during the compilation of the most re- 
liable atomic data. 

The abundances were varied until a satisfactory fit, as 
judged by visual inspection, was achieved for each element. 
We then changed the abundances in several steps of 0.1 - 0.2 
dex, until the fit was clearly worse. The resulting difference 
was used as a conservative estimate for the abundance error, 
or for an upper limit if no line was identified. Table[2]lists the 
lines used in this procedure, although not all lines could be 
used for all four stars. The best-fit models to the COS obser- 
vations are illustrated in Figs. [3] and [H and the metal abun- 
dances of the four white dwarfs are given in Table[3] (along 
with the previous a bundance studies were carried out fo r 
GALEX 1931+0117. (Vennes et al.|[2011bl : lMelis et al.ll201ll l. 
Notably, upper limits for N were always larger than solar 
relative to C. For Na, Ti, V, Mn (and additionally Ca in 
PG 0843-^516 as weU as Ca, Al, P, S, Ni in PG 1015+161) 
the upper limits were larger than solar relative to Si. We 
have used these (solar) values in the models, but it did not 
change the atmosphere structure and the results for the de- 
tected elements. 



3.2.1 Interstellar line absorption and airglow 

In all objects interstellar absorption is visible in the reso- 
nance lines of C ii, N i, O i. Si ii, and S ii. In SDSS 1228+1040, 
PG 1015+161, and GALEX 1931+0117 the interstellar ab- 
sorption lines are shifted blue-wards with respect to the 
photospheric lines by velocities of i; = 57, 36, and 
61 kms-\ respectively. In PG 0843+516, \v\ < 7 kms'^ 
and the interstellar lines are not fully separated from 
the photospheric features. However, the presence of some 
interstellar absorption is obvious from the line ratio of 
Cll 1334.5 A/C II 1335.7 A (Fig. |3]&|4}. Because the lat- 
ter line originates from a level only 0.008 eV above the real 
ground state, it is equally populated in a stellar photosphere, 
but not in the interstellar medium, where the blue compo- 
nent is much stronger in spite of a lower transition probabil- 
ity. Nevertheless, the abundances of C, O, Si, and S are ro- 
bust, as a sufficient number of excited transitions are present 
in the photospheric spectrum (Table[2]). 

The COS pipeline does not correct for airglow emission. 
Therefore, the reduced COS spectra can contain geocoronal 
lines of 1 1302, 1305, 1306 A whose intensity, and, to a 
lesser extent, profile shape, vary as a function of HST's or- 
bital day/night, and weakly with the Earth-limb angle. Air- 
glow is clearly seen in the spectrum of GALEX 1931+0117 
(Fig. 131 right panel), which affects the fit to the photospheric 
O I and Sill lines in this region. For Si, this is a minor prob- 
lem as there are many additional lines of Sill-iv. For O, 
another strong line in the COS spectra is Ol 1152 A. 



3.2.2 Silicon 

We notice relatively large differences of the silicon abun- 
dance determined from optical versus ultraviolet spectra in 



SDSS 1228+1040 and GALEX 1931+0117, for the latter also 
the oxygen abundances show this difference. There are at 
least three possible explanations: 

Uncertain atomic data. This is a perennial prob- 
lem, as there are many, and large differences in various 
compilations of atomic data. The O I resonance lines in 
GALEX 1931+0117 are perturbed by airglow, interstellar 
absorption and overlapping Sill lines (see above), and the 
ultraviolet abundance determination rests largely on one ex- 
cited line at 1152.1 A. Similarly, the optical O abundance 
is measured only from the 1 7777 A triplet jVennes et al.l 
I2OIOI : iMelis" et al.|[201lh . However, our abundance measure- 
ments for Si us e many lines in the u ltraviolet. In the recent 
compilation bv IBautista et al] l|2009l ) the authors combined 
several different computational methods, previous theoret- 
ical calculations by other authors, and experimental data 
into a "recommended" value for log gf. These values agree 
fairly well with the ultraviolet data from NIST that we have 
used. However, for the five optical lines they consider, the 
values are 0.25 — 0.30 dex smaller, though with errors as large 
as 0.3 dex. Using these values would increase the abundance 
determined from optical spectra, contrary to what would be 
needed for a more consistent solution. In addition, in a re- 
cent analys i s of u ltraviolet spectra for the DBZ star GD40, 
lJura et all (|2012l ) find a discrepancy between optical and 
ultraviolet abundances for Si of the same size, but in oppo- 
site direction - the abundances are smaller for the optical 
determinations. Since that study used the same models and 
atomic data as the one presented here, there is no indication 
that the atomic data are behind this discrepancy. 

Abundance stratification. Contrary to DB stars like 
CD 40 at similar temperatures, there are no convection zones 
in the atmospheres and envelopes of our four objects, which 
would act as a homogeneously mixed reservoir in the accre- 
tion/diffusion scenario. Assuming a steady state between the 
two processes, we thus expect a stratified abundance config- 
uration. Whether this can explain the observations will be 
studied in Sect. l3.2!3l 

Genuine variation of the accretion rates. As will also 
be discussed in the next section, the time scales for dif- 
fusion in these atmospheres are of the order of days. If 
the accretion rate is not constant the observed abundances 
may change on the same short time scales. Given that the 
COS and ground-based observations that we analysed were 
taken months to years apart, such variations can not be ex- 
cluded. Noticeable variations of the Call equivalent widths 
in the debris disc white dwar f G29-38 were reported by 
Ivon Hippel fc T hompso nI ll2007f). However, a simil ar study 
on the same star by Dc bes fc Lopez-MoralesI l|2008l ) did not 
find any variations in the line strengths. Thus, the cur- 
rent evidence for accretion rate variations on time scales 
of months to years is ambiguous, and a second-epoch COS 
observations of the stars studied here would be desirable. 

We also noticed an unidentified absorption feature be- 
tween 1400 and 1410 A, with a strength roughly correlated 
with the Si abundances. Such a feature has been discussed 
in the literature and related to an autoionisation line of Si 11 
or to a resonance feature in the p hotoionisation cross sec- 
tion jArtru fc Lanzll 19871 : iLanz et al. 1996, ). We have tested 
such a hypothetical line with their data for the oscillator 
strength and line width data. However, the width (« 80 A) 
is much too broad to lead to visible features in the spectrum. 
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Table 3. Element number abundances log[Z/H] and limits determined from the analysis of our HST /COS and optical spectra. For Si, 
we separately report the abundances determined from the optical data alone ("opt", see Sect. l3.2T2ll . and for Ca and Mg, we report both 
the abundances from a homogenous and a stratified atmosphere ("strat", see Sect. l3.2~3) l. the latter ones are preferred as t he stro ngest 
lines for bo t h elem ents are observed are in the optical. For GALEX 1931+0117, we also list the results of IVennes et al.l l l2011bl ) and 
iMelis et al.l l l201lf) . 



Element 


PG 0843+516 


PG 1015+161 


SDSS 1228+1040 
this paper 




GALEX 1931+0117 
Vennes et al. 


Melis et al. 


C 


-7.30 ± 0.30 


< -8.00 


-7.50 + 0.20 


-6.80 + 0.30 


< -4.15 


< -4.85 


O 


-5.00 + 0.30 


-5.50 + 0.20 


-4.55 + 0.20 


-4.10 + 0.30 


-3.62 + 0.05 - 


■3.68 + 0.10 


Mg 


-4.90 + 0.20 


-5.30 + 0.20 


-5.10 + 0.20 




-4.42 + 0.06 - 


-4.10 + 0.10 


Mg (strat) 


-5.00 + 0.20 


-5.30 + 0.20 


-5.20 + 0.20 








Al 


-6.50 + 0.20 




-5.75 + 0.20 


-6.20 + 0.20 






Si 


-5.20 + 0.20 


-6.40 + 0.20 


-5.20 + 0.20 


-4.75 + 0.20 






Si (opt) 






-4.70 + 0.20 




-4.24 + 0.07 - 


■4.35 + 0.11 


P 


-6.60 + 0.20 




< -7.30 


-7.00 + 0.30 






S 


-5.50 + 0.30 




< -6.20 


-6.60 + 0.20 






Ca 




-6.30 + 0.20 


-5.70 + 0.20 




-6.11 + 0.04 - 


-5.83 + 0.10 


Ca (strat) 




-6.45 + 0.20 


-5.94 + 0.20 








Cr 


-5.80 + 0.30 


< -5.80 


< -6.00 


-6.10 + 0.30 




■5.92 + 0.14 


Mn 












-6.26 + 0.15 


Fe 


-4.60 + 0.20 


-5.50 + 0.30 


-5.20 + 0.30 


-4.50 + 0.30 


-4.43 + 0.09 - 


-4.10 + 0.10 


Ni 


-6.30 + 0.30 




< -6.50 


-6.70 + 0.30 




< -5.60 



We have also included the Si ll photoionisation cr oss sec- 
tions from the Opacity Project (ISeaton et al.|[l993 ). which 
indeed show a resonance maximum in this spectral region. 
But again, the Si abundance is too small to let this feature 
show up in the spectrum. 

Our model uses the six Sin lines at 1403.8, 1404.2, 
1404.5, 1409.1, 1409.9, and 1410.2 A in this range (Table[2l). 
The first two have the source "guess" in VALD, the first 
three have no entry in NIST, and the log gf values of the 
strongest line (1410.2 A) differ by ~ 0.8 dex between the two 
databases. The upper levels of the transitions have a par- 
ent configuration belonging to the second ionisation limit of 
Sill. They are still ~ 0.7 eV below the first ionisation limit 
and thus not strictly auto-ionising. However, the broaden- 
ing may well be underestimated by our simple approxima- 
tion formulae. In summary, the atomic data of the lines in 
the region are very uncertain and may be the explanation 
for the broad feature. However, with the present data we 
cannot prove that hypothesis. 

Finally, we note that the Si iv 1394,1403 A doublet in 
PG 0843+516 is very poorly fit by our atmosphere model 
(Fig. El- A weaker additional Siiv 1394,1403 A absorption is 
also seen in the spectrum of SDSS 1228+1040 (Fig.lHl. We 
interpret this as evidence for absorption by hot gas close to 
the white dwarf, see the discussion in Sect.[71 



3.2.3 Diffusion and stratified atmosphere models 

In the absence of a convection zone there is no deep homoge- 
nous reservoir in our DAZ sample, and therefore there is no 
straightforward definition of diffusion time scales. Adopting 
the usual definition, i.e. dividing the mass of some element 
above a layer in the envelope or atmosphere of the star by the 
diffusion flux, results in diff usion time scales that str ongly 
depend on the chosen layer. iKoester fc WilkenI l|2006l ) and 
iKoestej (,2009 ) defined the Rosseland optical depth t = 5 as 



the "standard" layer, assuming that no trace of any heavy 
element below this would be seen in a spectrum. 

However, a more consistent way to determine the abun- 
dances in the accreted material, which is the quantity ulti- 
mately desired, is the assumption of a steady state between 
accretion and diffusion throughout the whole atmosphere. 
At Rosseland optical depth r = 2/3, and typical conditions 
for the observed ultraviolet spectra, the diffusion times in 
the four white dwarfs analysed here are ~ 0.4 to four days. 
Assuming that the accretion rate does not vary over such 
time scales, we can use the condition of constant fiow of an 
element with mass fraction X{t) 

pXv = const (1) 

with p and v the mass density and the diffusion velocity 
of this element, p and v are known from the atmosphere 
model and diffusion calculations, and X{t = 2/3) is derived 
from the spectral analysis. This determines the diffusion flux 
at T = 2/3. In steady state, as it is the case for the DAZ 
analysed here, the diffusion flux is constant throughout the 
atmosphere, and is equal to the accretion rate polluting the 
atmosphere. The constant diffusion flux then in turn allows 
the determination of the abundance stratification X{t) (see 
also IVennes et al]|2011bl . for a thorough discussion). 

We calculated new stratified models and synthetic spec- 
tra for all objects, using the steady state condition and the 
abundances (at t = 2/3) from Tabled The resulting spec- 
tra are almost indistinguishable from those of the homoge- 
neous atmospheres; the only exception are small increases of 
the optical Mgll and Call line strengths. The small change 
can easily be explained by the structure of the stratified 
atmosphere. In these models pv increases with depth, and 
consequently the abundance decreases. On the other hand a 
monochromatic optical depth of ~ 2/3 is reached in the ul- 
traviolet near Rosseland optical depth of tross — 2/3, while 
it is reached at tross ~ 0.15 for A = 4480 A, i.e. higher 
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Table 4. Diffusion fluxes pXv[gs^^] within the white dwarf atmospheres, which are equal to the rates at which planetary debris 
material is accreted. S gives the sum of the accretions rates of all detected ele ments (i.e . not in cluding those with upper li mits). For 
GALE X 1931+0117, we adopt in our calculation the Mg and Ca abundances of iVennes et al.l ll2011bl) and the Mn abundance of iMelis et al.l 



Element PG 0843+516 PG 1015+161 SDSS 1228+1040 GALEX 1931+0117 



c 


1.66 


X 


10'> 


4.65 


X 


10-* 


1.25 


X 


10^ 


4.57 X 10^ 


o 


9.27 


X 


10^ 


3.78 


X 


10^ 


2.70 


X 


10* 


5.61 X 10* 


Mg 


4.47 


X 


10^ 


2.66 


X 


10^ 


3.21 


X 


10^ 


1.47 X 10* 


Al 


2.09 


X 


10'5 








1.18 


X 


lO'^ 


3.08 X 10^ 


Si 


4.77 


X 


10^ 


3.64 


X 


10'' 


4.80 


X 


10^ 


9.93 X 10^ 


P 


2.44 


X 


10^ 








< 5.24 


X 


iqs 


7.57 X 10^ 


s 


3.92 


X 


10^ 








< 9.46 


X 


10^ 


2.64 X 10"^ 


Ca 








4.84 


X 


IQf* 


1.57 


X 


10^ 


8.10 X 10^ 


Or 


3.81 


X 


10^ 


< 3.85 


X 


10^ 


< 2.29 


X 


lO'^ 


1.37 X 10^ 


Mn 




















1.06 X 10^ 


Fe 


7.11 


X 


10* 


9.50 


X 


10^ 


1.72 


X 


iqs 


6.45 X 10* 


Ni 


1.66 


X 


107 








< 9.98 


X 


10^ 


4.71 X 10'3 


S 


1.02 


X 


109 


1.68 


X 


iqs 


5.61 


X 


iqs 


1.50 X 10^ 



in the atmosphere, where the abundance is correspondingly 
higher. 

For PG 0843+516, PG 1015+161, and SDSS 1228+1040, 
the Ca and Mg abundances were obtained from the optical 
data (Sect.[2I2J and our models. We have iterated them by 
fitting to stratified models (denoted with "strat" in Table[3]). 
For GALEX 1931 +0117, we adopted the photospheric Mg an 
Ca abun dances of [ Vermes et aU l|2011bl ) and the Mn abun- 
dance of lMelis et al.l ( 2011 ) to calculate the corresponding 
diffusion fluxes. 

As a result we have to conclude that diffusion and a 
stratified abundance structure lead only to minor adjust- 
ments of the abundances that cannot explain the large dis- 
crepancy between optical and ultraviolet determinations for 
silicon. There is, however, an important caveat to this con- 
clusion. Our diffusion calculations use only the surface grav- 
ity (and as a minor effect the t emperature gradient for 
thermal diffusion) as driving force. IChaver fc Dupuii (|2010l ) 
have recently demonstrated that for silicon, radiative levi- 
tation can lead to a negative effective gravity and support 
the atoms in the outer layers of the atmosphere against dif- 
fusion. They only published detailed data for a DAZ model 
with 20 000 K and \ogg — 8.00, and in their model only 
abundances smaller than log [Si/H] = —8.0 are really sup- 
ported, because the lines saturate at higher abundances, ef- 
fectively reducing the radiative support. However, it is quite 
feasible that even if the atoms are not totally supported, the 
diffusion velocity would be smaller, changing the abundance 
gradient. The answer to this puzzle will have to await sim- 
ilar, detailed models for a variety of stellar parameters and 
heavy elements that can be tested against the large range 
of Si abundances found in our snapshot survey (Gansicke et 
al. in prep). 

Other points worth mentioning are that the determi- 
nation of an effective io n charge with the sim ple pressure 
ionisation description of iPaguette et al.l (| 19861 ) is not ap- 
propriate in the absence of deep convection. We have used 
the usual Saha equation (with a small lowering of the ion- 
isation potential from non-ideal interactions) to determine 
the abundances of different ions from an element. The dif- 
fusion velocity is then calculated as a weighted average of 



the ionisation stages. This procedure was already used in 
iKoester fc Wfi kcn (2006) and Koestcr (2009) for the models 
without or with only a shallow convection zone, although not 
explicitly stated in those papers. New in our present calcu- 
lation is the consideration of neutra l particles, fol l owing t he 
discussion and methods outlined in lVennes et al.l l|2011ah . 

The main results of our calculations are the diffusion 
fiuxes, Xpv, for each element, which are assumed (in steady 
state) to be the abundances of the accreted matter. These 
are summarised for the four objects in Table 2] The total 
diffusion fluxes ( = accretion rates) are obtained by multi- 
plying t hese fluxes wi th iivR^^, where we used the cooling 
tracks of lWoodI ^99^) to obtain the white dwarf radii from 
Tcff and log g. The mass fluxes ( = accretion rates) of the in- 
dividual elements, as well as their sum, are shown in Fig. [5] 
and discussed in Sect. |S] The number abundances of the 
circumstellar debris are then calculated from the diffusion 
fluxes via 



N(X) _ M(X) A(Si) 



N(Si) M(Si) A(X) 



(2) 



where A is the atomic mass. The implications that these 
abundances have on the nature and origin of the circumstel- 
lar debris are discussed in detail in Sect .[5] 



4 NOTES ON INDIVIDUAL WHITE DWARFS 

In the following sections, we give a brief overview of previous 
work on the four white dwarfs that we have analysed, as well 
as a summary of the key results of our observations. 



4.1 PG 0843+516 

PG 0843+516 was identified as a DA white dwarf 
in the Palomar-Gr een Survey ijGreen et all Il986l ). and 
iLiebert et aU (|2005l ) obtained = 23 870 ± 392 K, log g = 
7.90 ± 0.05 from the analysis of a high-quality optical spec- 
trum. The best fit to our HST data was Tcff = 23 095 ± 
230 K, logg = 8.17 ± 0.06. Our COS spectrum reveals 
PG 0843+516 to be an extremely polluted DAZ white dwarf 



© 2006 RAS, MNRAS 000.[llfT6l 



10 B. T. Gdnsicke et al. 



r 




1* 














1 1 




r 




t 

1 










ft (( 1 


\ 




t : 














I 






( 










1 




1 




( 1 


r 








1 












• GALEX 1931+01 17 


















• PG 0843 + 516 




















• PG 1015+161 








1 




1 






1 


• SDSS 1228+1040 

1 1 1 1 1 


1 




c 





s 


Mn 


P 


Cr 


Si Mg Fe Ni 


Ca 


Ti Al 


Sum 



Figure 5. Accretion rates of tiie elements detected in our four 
targets. Tlieir sum is given in the riglit-most column. 



(Fig.[T] & El, with an accretion rate of ~ lO'gs ^, placing 
it head-to-head with GALEX 1931-^0117 (Sect.|6|. We iden- 
tified in the COS spectrum photospheric absorption lines 
of C, O, Al, Si, P, S, Fe, Cr, and Ni, plus Mg in the op- 
tical WHT spectrum, the second largest set of elements 
detected in a DAZ white dwarf. The fact that the metal 
pollution of PG 0843-1-516 went unnoticed in the published 
high-quality intermediate resolution spectroscopy underlines 
the strength of our ultraviolet survey for young and rel- 
atively w arm white d warfs accreting planetary debris. We 
note that IXu fc Jural |2012il recently detected infrared flux 
excess at PG 0843-1-516 in an analysis of archival Spitzer 
data, making this the second white dwarf (af ter G29-38, 
IZuckerman fc Beck"iinlll987l : iKoester et al.lll997l ') where cir- 
cumstellar dust was found without prior knowledge of pho- 
tospheric metal pollution. 



4.2 PG 1015+161 

PG 1015-1-161 is another DA white dwarf d is covered in the 
Palom ar-Green Survey jCreen et al.l 1 19861 ). iLiebert et al.l 
l|2005l ) determined T^s = 19 540 ± 305 K, log g = 8.04 ± 0.05 
from optical spectroscopy. High-resolution spectroscopy of 
PG 1015-1-161 was obta ined as part of the SPY p roject 
iNapiwotzki et all [ioOlI 'l. from which iKoester et all (|2009l ') 
measured Tea = 19 948 ± 33K and logp = 7.925 ± 0.006. 
Our fit to the HST spectrum give s in T^ s = 19 200 ± 180 K, 
log g = 8.22 ± 0.06. IKoester et akl (120051 ) detected of a pho- 
tospheric CallK absorption line in the SPY data, with 
a number abundance log [Ca/H] — —6.3, which triggered 
follow-up observations wit h Spitzer that r evealed the pres- 
ence of circumstellar dust (|Jura et al.ll2007i ). The COS spec- 
trum contains absorption lines of O, Si, and Fe. In addition 
to CallK, we detected Mgll 4482 A in the SPY spectrum. 
PG 1015-1-161 has the lowest accretion rate among the four 
stars discussed in this paper. 



4.3 SDSS 1228+1040 

lEisenstein et al.l (|2006l ) identified this DA white dwarf in 
Data Release 4 of the Sloan Digital Sky Survey, and found 



Tcff = 22 125 ± 136 K, logq = 8.22 ± 0.02 from a fit 
to the SDSS spectrum. iGansicke et al.l l|2006l ) discovered 
double-peaked emission lines of Can 8498,8542,8662 A as 
well as weak Fe ll emission lines and Mg ll 4482 A absorp- 
tion, and concluded that SDSS 1228+1040 accretes from a 
volatile-depleted gaseous circumstellar disc. The Call lines 
form in a region extending in radius from a few tenths 
R© to ~ 1.2 R0, no emission is detected from closer in 
to the white dwarf (but see Sect. [7|). Spitzer observations 
show ed that SDSS 1228+104 also exhibits an infrared ex- 
cess (|Brinkworth et al.|[2009l ). and that there is a large ra- 
dial overlap between the gaseous and dusty components of 
the disc. Yet, the strong Can emission li nes require a gas 
temp erature of T ~ 4000 - 6000 K (e.g. iHartmann et all 
l201ll ). substantially exceeding the sublimation temperature 
of the dust. This implies the thermal decoupling of the 
gas and dust, most likely in the form of a complex verti- 
cal temperature structure, with hotter, optically thi n gas 
on top cooler, pr obably optically thick dust (|Kinneaill201ll : 
iMelis et al.ll2010l ). Irradiation from the whit e dwarf is suf- 
ficient to explain this temperature inversion l|Kinneaill20"Tll : 
iMelis et al.ll201 0'). but the origin of the gas found at radii 
larger than the sublimation radius is uncle ar, and may be 
related to relatively fre sh disruption events (|Gansicke et al.l 
I2OO8I : iMelis et al]|2010l) or the intrins i c evolution of the de - 
bris disc (|Bochkarev fc Rafi"kwll201ll : iMetzger et al.ll2012l ). 
Among the four white dwarfs studied here, SDSS 1228+1040 
is the only one that exhibits emission lines from a gaseous 
disc. 

The COS spectrum of SDSS 1228+1040 contains ab- 
sorption lines of C, O, Al, Si, Cr, and Ni. SDSS 1228+1040 
was observed outside the snapshot program described in 
Sect. 12. ll and our COS spectroscopy extends up to 1790 A, 
i.e. 360 A further than that obtained for the other three 
white dwarfs. This extended wavelength range includes ad- 
ditional strong lines of Sin, Alll, and Aim, but no further 
elements. Our high-quality average UVES spectrum is used 
to determine the abundances of Mg and Ca, bringing the 
total number of detected elements in SDSS 1228+1040 to 
eight. 

We fitted the SDSS spectrum, finding T^s = 22 410 ± 
175 K, log (7 — 8.12 ± 0.03, whereas a fit to the ultraviolet 
spectrum gives Toff = 20 565 ± 82 K, logg = 8.19 ± 0.03. 
This discrepancy underlines that, for high-quality data, 
the uncertainties are dominated by systematic rather than 
statistical errors. As a compromise we take the weighted 
mean of the latter two results with increased errors, T^g = 
20 900 + 900 K, logg = 8.15 + 0.04. 



4.4 GALEX 1931+0117 

As part of a spectro scopic identification program of 
ultraviolet-excess obiects lVennes et al.l l|201(]| ) recently iden- 
tified GALEX 1931+0117 as a nearby (~ 55 pc) D AZ white 
dwarf IVennes et all (|2010l ) and lMelis et all (|201ll ) analysed 
optical spectroscopy, and obtained Tefi = 20 890 ± 120 K, 
logg = 7.90+o;°g and Tcft = 23 470 + 300 K, logg = 
7.99 + 0.05, respectively. Our best-fit parameters from the 
HST /COS spectrum are ns = 21 200 ± 50 K, l o gg = 
7.91 ± 0.02, consistent with that of IVennes et al] (|2010l ) 
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but somewhat lower than that of lMeHs et al. (12011 The 
VLT/U VES spectroscopy obtained bv lVennes et'aLri|2010l. 
l2011bl ) revealed strong metal lin es of O, Mr, Si, Ca , and 
Fe, indicating ongoing accretion. IVennes et all (|2010l ) also 
showed that the 2MASS H- and A'-band fluxes exceeded 
those expected from the white dwarf, and suggested a close 
brown dwarf or a dusty d e bris d isc as origin of the accret- 
ing material. iDebes et al.l l|201ll ) ruled out the presence of 
a sub-stellar companion based on the infrared fluxes de- 
tected by WISE, and argued that the white dwarf accretes 
from a dusty disc. This was independently c onfirmed by 
VLT/ ISAAC near-IR observations obtained bv lMelis et al.l 
l|201ll ). who also measured abundances for Cr and Mn. 

Our HST /COS spectroscopy provides independent 
measurements for O, Si, Cr, and Fe, as well as the first de- 
tection of C, Al, P, S, and Ni, bringing the total number of 
elements observed in the photosphere of GALEX 1931-1-0117 
to 11 (TableO. As discussed in Sect.[331 the O, Si, Cr, and 
Fe abundances that we derive from the COS spec t roscopy 
are lower than t hose determined by IVennes et al.l l|2011bl ) 
and iMelis et al.l (|201ll ). However, the discussion of the na- 
ture of the planetary material is us ually based on rela - 
tive metal-to-metal abundance ratios (|Nittler et al.|[20o3 ). 
which are more robust than absolute abundances measure- 
ments. Figure[6] compares the metal abundances determined 
for GALEX 1931-1-0117 normalised with respect to Si, and 
relative to the corresponding ratios for the chemical compo- 
sition of the bulk Earth. It is evi dent that ou r meta l-to-Si ra- 
tios are consistent with those o f lMelis et all l^2011^■ whereas 
the Mg/Si, Fe/Si, and Ca/Si ratios of IVennes et all l|2011bl ') 
are systematically lower. 



5 THE NATURE AND ORIGIN OF THE 
CIRCUMSTELLAR MATERIAL 

The four white dwarfs studied here have diffusion time scales 
of a few days (Sect.[3I23Ji and we can therefore safely as- 
sume that we observe them in accretion-diffusion equilib- 
rium. In other words, the abundances of the circumstellar 
debris can be determined from the photospheric analysis 
without any additional assumptions regarding the history 
of the accretion rate that are n ecessary for stars with very 
long diffusion time scales (e.g. iKlein et al.ll201ll ). In what 
follows, we discuss the abundances of the circumstellar de- 
bris normalised to Si, the main rock- for ming element, as is 
common use for solar-system objects (e.g. lLodders fc FeglevI 

Figure[6] (right panel) illustrates the metal-to-Si ratios 
of the planetary debris around the four white dwarfs rela- 
tive to the same a bundances of the bulk Earth model by 
iMcDonouglil (|2000l ). The first striking observation is that 
the C/Si ratios of all four stars (including one upper limit) 
are much lower than that of CI chondrites, and in fact agree 
within their errors with the C/Si value of the bulk Earth 



^ iMelis et al.l ll201lh discuss the discrepancy between their model 
and the GALEX fluxes. From their Tablet, it appears that 
they did not correct for the non-linearity of the GALEX detec- 
tor s for bright targets. The corrected GALEX magnitudes given 
bv lVennes et al] ll201Cll ) are in good agreement with our best-fit 
model. 



model. While the C abundance of the bulk Earth is subject 
to some model-dependent assumptions (see the left panel 
of Fig.E] for an alte rnative chemical model of the Earth by 
lAllegre et alll200ll ). these uncertainties are comparable to 
the errors in our abundance determinations. 

For comparison, we include in Fig. |6] the abundance ra- 
tios of three white dwarfs that accrete from the wind of a 
close M-dwarf companion, that were also observed as part of 
our COS snapshot programmfl. The only elements detected 
in the COS spectra of these three stars are C, O, Si and S, 
and they exhibit high abundances in C and S, as expected 
for the accretion of solar-like material. The extremely low 
abundances of the volatile C found for the debris around 
the four white dwarfs strongly underlines its rocky na ture. 
This corro b orate s the previous studies of lJural (I2006D and 
lJura et all l|2012l ). who found strong evidence for substan- 
tial depletion of C around three DB white dwarfs. 

However, Fig. [6] also shows that there is a significant 
scatter among the individual abundances for a given ele- 
ment. Among the four targets, the abundances of the debris 
in SDSS 1228+1040 most closely resembles those of the bulk 
Earth. PG 10154-161 stands out by having all detected ele- 
ments over-abundant with respect to Si, when compared to 
the bulk Earth. An interesting trend is seen in PC 0843-1-516, 
where Fe, Ni, and S are significantly over- abundant, and, in 
fact, broadly consistent with the abundance ratios of the 
core Earth model. In particular, the volatile S is extremely 
overabundant with respect to C, compared to the bulk sili- 
cate Earth. In melts, S will form FeS, and hence be depleted 
from remaining minerals. The affinity of S to Fe is thought 
to be the reason for the depletion of S in the silicate mantle 
of the Earth, as it will ha ve settled int o the Earth's core 
in the form of iron sulfide l|Ahrenslll979l : iDreibus fc Palmd 
11996). Similarly, also Cr is significantly over-abundant in 
PG 0843-1-516 with respect to the bulk Earth. While Cr is 
a moderately volatile element, the depletion of Cr in the 
silicate Earth is thought to be due to partitioning into the 
Earth's core l|Movnier et al.l [20111 ). Finally, the refractory 
lithophile Al is under-abundant compared to the silicate 
Earth. Thus, the abundance pattern seen in PG 0843+516 
suggests that the planetary debris is rich in material that 
has undergone at least partial melting, and possibly differ- 
entiation. A possible test of this hypothesis would be a mea- 
surement of the abundance of Zn, a lithophile element with 
a simi lar volatility as S that is not depleted into iron melt 
(|Lodde ra 2003). and it will be important to test whether the 
refractory lithophile Ca is depleted at a similar level as Al. 
The most promising feature to measure the Zn abundances is 
the Zn II 2026,2062 A resonance doublet, and CallK should 
be easily detectable in high-resolution optical spectroscopy. 

To further explore the chemical diversity of the plan- 
etary debris around the four white dwarfs studied here, 
we compare pairwise a range of metal-to-Si abundance ra- 
tios with those of the bulk Earth and bulk silicate Earth 
dMcDonou gh 2000'), as well as with tho se of a variety of me- 
teorites (taken from lNittler et al]|2004f ). We inspect first the 



A more detailed discussion of these binaries will be pubUshed 
elsewhere. Here, they merely serve as "abundance standard white 
dwarfs" which accrete material with abundance ratios that are 
expected to be close to solar, i.e. rich in volatiles. 
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Figure 6. Heavy clement abundances derived for th e circumstellar deb ris at the four white dwarf targets (see Table |4]l, relative to Si, 
and normalised to the same ratios of the bulk Earth llMcDonougbll200ClV The elements are arranged, left to right, in order of increasing 
sublimation temp erature I ILodd crs 2003) . Left p anel: abundances for the debris around GALEX 1931+0117 (green: this paper, blue: 
iMelis et al.ll20l"ll . maeenta: IVennes et al]|2011bl) . Also shown are the abundance ratios for the core Earth (which makes up ~ 1/3 of 
the Earth's mass, open triangle s) and the silicate Earth (i.e crust and mantle, which make up ~ 2/3 of the Earth's mass, open circles). 
The bulk Earth composition of lAUegre et al.| fl995) is shown as solid black squares, illustrating the level of uncertainty in the (model- 
dependent) composition of the Earth. The sho rt-dashed line shows the abundance ratios of CI chondrites, the long-dashed line those 
corresponding to solar abundances (both from iLodderj [ioOSi) . Right panel: Metal-to-Si ratios for PG0843-I-516 (blue), PC 1015-1-161 
(magenta), SDSS 1228-1-1040 (red), and GALEX 1931-1-0117 (green). Shown in orange are the abundance ratios of three white dwarfs 
that accrete from the wind of a close M-dwarf companion. As expected for the accretion of material with near-solar abundances, the 
volatiles C and S are found to be strongly enhanced compared to the four white dwarfs that host debris discs of exo-terrestrial material. 



relative abundances of Al and Ca, which are two of the three 
most abundant refractory lithophile elements (the third one 
being Ti), i.e. elements that sublimate only at very high 
temperatures, and that do not enter the core in the case of 
differentiation. Therefore, the Al/Ca ratio is nearly constant 
across most classes of meteorites, and hence, the Al/Si values 
determined from many solar-system bodies follows a linear 
correlation with Ca/Si (Fig.[7l top right). Finding that the 
abundances for the debris discs, where Al, Ca, and Si are 
available, generally follow that trend is reassuring, as large 
variations in the relative Al and Ca abundances would cast 
doubts on the overall methodology using white dwarf pho- 
tospheres as proxies for the abundances of the circumstellar 
material. 

The relative abundances of O, Si, Mg, and Fe, which 
are the major constituents of the terrestrial planets in the 
solar system, show substantial variations between different 
meteorite groups (Fig. [71 top left and bottom right panels), 
and at least as much scatter between the individual white 
dwarfs. The difficulty with these elements is that they form a 
range of different minerals (metal oxides) , depending on the 
prevailing pressure and temperature. Iron in particular may 
occur as pure metal, alloy, or mineral, and is subject to dif- 
ferentiation into planetary cores. Oxygen, on the other hand, 
can be be locked in a wi de range of oxides (see the discus- 
sion by Klein et al. 2010l), or potentially water (iKlein et al.l 
I2OIOI : lJura fc XulbOldlFarihi et al.ll201ll : I Jura fc Xull2012h . 
Therefore, the relative abundances of O, Si, Mg, and Fe will 
vary according to the processing that material underwent 
(e.g. condensation, melting, and differentiation), and it is 
maybe not too surprising to find that the debris around 
white dwarfs exhibits at substantial degree of diversity, as 
it represents different planetary systems formed around dif- 



ferent stars. We note that the debris at PC 08434-516 falls 
close to the abundance ratios of Pallasites, a class of stony- 
iron meteorites. This further supports our hypothesis that 
PG 0843-1-516 is accreting material in which iron has under- 
gone (partial) melting. 

Another interesting pair of elements is C and O (Fig. [71 
lower left panel). The possible range of the C/O ratio 
among exo-planets has been subject to intense discussion. 
It is thought that for C/O > 0.8 in the proto-planetary 
discs, the ambient chemistry will favour solid "carbon plan- 
ets", that are dominated by carbides rather than oxides 
llKuchner fc Se aacr 2005). The possible existence of carbon 
planets has gained some support by the recent report of a 
C/O value exceeding unity i n the atmosphere of the tr ansit- 
ing hot Jupiter WASP-12b ([Madhusudhan et al.l^2011^ . and 
by abundance studies that found a significant fraction of 
exo-planet host stars having C/ O > 0.8 (Petigura & Marcvl 
l201ll : lDelgado Mena et "alll2010l ). but see lFortnevI (|2012l ) for 
a critical discussion. 

Planetary debris at white dwarfs provides a unique op- 
portunity to probe the C/O ratio of exo-terrestrial ma- 
terial. However, measuring C abundances in white dwarfs 
is challenging, as the optical detection of carbon in cool 
white dwarfs is usually related to d redge-up from the 
core rather than external pollution (e.g. [Pufour et al.ll2005l : 
iKoester fc KnistI I2OO6I : [Pesharnais et all |2008| ) . At higher 
temperatures, where convective dredge-up can be excluded, 
suitable lines of C are only found at ultraviolet wavelengths. 
As mentioned above, the four stars studied here have very 
similar (low) C/Si ratios, but do show a range of O/Si ra- 
tios. Nevertheless, the debr is around all fou r stars stud- 
ies here, as well as CD 40 (|jura et al.l [201 have —3 < 
log(C/0) < —2.3, very similar to the bulk silicate Earth, 
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Figure 7. The chemical abundances of planetary debris, determined from the photospheric studies of polluted white dwarfs, reveal a large 
degree of diversity. The four panels illustrate a range of metal-to-Si n umber abun dance ratios for the four stars analysed in this paper, 
compared to t hose of bulk Earth and bulk silicate Earth (BE and BSE: iMcDonoughi . 2000). solar abundances and CI chondrites (S and CI; 
iLodder j|2003l ). and several meteorite clas ses (gray = carbona ceous Chondrites, green = Mesoderites, blue = Pallasites, red = Diogenites, 
orange = Howardites, magenta = Eucrites: [ Nittler et al]|2004l). Also sho wn, in light blue, are the abundance r atios for the polluted DB 
white dwarfs GD362 llZuckerman et al.ll2007t ). CD 40 jjura et alJlioT^) . and HS 2253+8023 l lKlein et al.ll201lh . 



log(C/0) ~ —2.5, and are hence representative of solar sys- 
tem minerals. 



6 ACCRETION RATES 

Estimating accretion rates for metal-polluted white dwarfs 
is notoriously difficult, as it is based on scaling from the 
elements detected in the photosphere to an assumed bulk 
composition of the accreted material. In addition, in the case 
of white dwarfs with significant convective envelope masses, 
only the average accretion rate over the diffusion time scale 

can be obtained. 

iKoester fc WilkenI l|200d ) calculated accretion rates for 
38 DAZ white dwarfs based on the abundance of Ca, 
and adopting solar abundances for the accreting mate- 
rial. For PG1015-fl61, these assumptions implied M ~ 
2 X lO'^^gs"^. Since then, it has become increasingly clear 
that many, if not most, metal-polluted (single) white dwarfs 
accrete vol atile-depleted mate rial from circumstellar plane- 
tary debris. iFarihi etal] l|2009l ) estimated accretion rates for 



53 metal-polluted white dwarfs following the prescription of 
iKoester &: WilkenI {200^) , but scaling the results by the typ- 
ical gas-to-dust ratio in the interstellar medium to account 
for the absence of H and He in the accreted debris, resulting 
in M ~ 2 X 10^ g s"^ for PG 1015-^161. 

The uncertainty in the estimated accretion rates can be 
greatly reduced if photospheric abundances for the major 
constituents of the debris material can be measured. While 
we do not detect all elements that are likely present in the 
circumstellar debris at the four white dwarfs studied here, 
we have determined the accretion rates of all the major el- 
ements, in particular O, Si, Mg, and Fe (Sect . 13^2^ . The 
accretion rates of all detected elements, as well as their 
sum are given in Table |4l and are illustrated in Fig. (5] 
For PG 1015-^161, we find M ~ 1.7 x lO^gs'S which is 
strictly speaking a lower limit, however, the undetected el- 
ements (e.g. Al, S, Ti, Mn, Cr) are unlikely to contribute 
more than 10% of the total accretion rate. Similarly, we find 
the accretion rates of PG 0843-^516, SDSS 1228-^1040, and 
GALEX 1931-^0117 to be M ~ 1.0x10^ gs'S 5.6x10** gs'S 
and 1.5 x 10^ gs"'^, respectively. 
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7 HOT CIRCUMSTELLAR GAS 

The discs around white dwarfs are passive, i.e. their emis- 
sion is solely due to the thermal reprocessing of inter- 
cepted stellar flux. The inner disc radius where typical 
dust grains will rapidly subli mate is determined by t he lu- 
minosity of the white dwarf (jvon Hippel et allbOOTl ). The 
gaseous material will viscously spread, both flowing in- 
wards onto the white dwarf, and outwards over the dusty 
disc, potentially accelerating the inwards migration of the 
dust via aerodynamic drag ()Rafikovll201ll ). While gaseous 
material orbiting at radii coincident with circumstellar 
dust is observed in a number of systems in the f o rm of 



2007 



2012 



doublc-pcakcd omission lines ("Gansicke et al.' '200' 
'2OO8; Brinkworth ot al. 2009, 2012; Mchs ct al. 2011] 
Farihi et all l2012l : iDufour et al.ll2012i ). there has yet been 
no detection of gaseous material well inside the sublimation 
radius. 

Inspection of Fig. [3] reveals that the strength of 
the Si IV 1394,1403 A doublet in PG 0843-1-516 is extremely 
under- predicted by the photospheric model. These Si iv lines 
correspond to the highest ionisation energy of all tran- 
sitions detected in the COS spectrum. For the tempera- 
ture and the Si abundance of PG 0843-1-516, the observed 
strength of the Siiv lines is absolutely incompatible with 
a purely photospheric origin. The most plausible explana- 
tion is that there is additional absorption along the line 
of sight, associated with hot gas close to the white dwarf 
that is optically thin except for the strong resonance lines 
of high-ionisation species, such as Silv. In fact, extremely 
similar features were found in the far-ultraviolet observa- 
tions of cataclysmic variables, i.e. white dwarfs that ac- 
crete from a (hydrogen-rich) accretion disc that is in turn 
fed by Roche-lobe overflow of a close M-dwarf companion. 
HST/GYmS and FUSE spectroscopy of the white dwarf in 
UGem contains very strong absorption of Nv 1239,1243 A 
and Ovi 1032, 1038 A that can not form in the ~ 30 000 K 
photosphere, as well as excess absorption in Si II 1394,1403 A 
i|Sion et al.lll998l : iLong fc Gillilandlll"999l : iLong et~ai]|2006h . 
All three high-ionisation doublets are red-shifted with re- 
spect to the systemic velocity of the white dwarf, but some- 
what less so than the lower-ionisation photospheric lines, 
which are subject to the gravitational redshift at the pho- 
tospheric radius. These observations were interpreted as ev- 
idence for a hot 80 000 K) layer of gas sufficiently close 
to the white dwarf to still experience a noticeable gravi- 
tational redshift. Measuring the central wavelengths of the 
strong Si iv 1394, 1403 A lines in PG 0843-1-516, we find that 
they are blue-shifted with respect to the photospheric fea- 
tures by ~ 25kms~^, which implies a height of ~ 1.5 white 
dwarf radii above the white dwarf surface. This assumes that 
there is no significant flow velocity, which seems reasonably 
well justified given the symmetric shape of the Si IV profiles. 

A discrepancy between the best-fit white dwarf model 
and the region around the Silv doublet is also seen in 
the COS spectrum of SDSS 1228-H1040 (Fig. H bottom left 
panel), however, in this star, the additional absorption is 
rather weak. These additional absorption features are clearly 
blue-shifted with respect to the photospheric lines, however, 
the relatively low signal-to-noise ratio of the spectrum pre- 
vents an accurate determination of this offset. 

For PG1015-M61 and GALEX 1931-^0117, the photo- 



spheric fits match the observed Silv lines well, i.e. there 
is no evidence for any additional absorption component. 
Given that these two stars have, respectively, the lowest 
and highest accretion rate of our small sample (Sect.|6]), 
there seems to be no clear correlation between the detec- 
tion of absorption from highly ionised gas to the mass flow 
rate onto the white dwarf. A key difference between the two 
stars where circumstellar Si iv absorption is detected is that 
SDSS 1228+1040 also shows strong emission lines from cir- 
cumstellar gas, which indicate a relatively high inclination of 
the accretion disc. In contrast, no gaseous emission is found 
in PG 0843-1-516 (Gansicke et al. in prep). Identifying addi- 
tional absorption features from these hot layers of gas would 
provide substantial constraints on the physical parameters 
in the corresponding regions. The strongest line seen in cat- 
aclysmic variables, N v, is naturally absent in the white 
dwarfs accreting rocky de brifl but the O vi 1032,1038 A 
doublet detected in U Gem (|Long et al.ll200d) is a promising 
candidate. 



8 CONCLUSIONS 

Recent years have seen a surge of interest in the evo- 
lution of extra-solar planetary systems th rough the late 



phase s in the lifes of their host st ars (e.g. Burleigh et al.l 



2002'; Dobcs & Sigurdsson '2002'; "Villaver & Livio' |2007| . 
2009; Nordhaus__et_al, 2010; Pi Stefano__ct_aL 2010). While 
no planet has yet been discovered or biting a white dwarf 
ijHogan et al.1 l2009l : iFaedi et al.l I2OIII ) , significant progress 
has been made in the discovery and understanding of plan- 
etary debris discs around white dwarfs. 

Our COS study substantially increases the number of 
polluted white dwarfs for which a wide range of chemical 
elements have been detected. We find that the C/Si ratio 
is consistent with that of the bulk Earth, which confirms 
the rocky nature of the debris at these white dwarfs, and 
their C/0 values are typical of minerals dominated by Fe 
and Mg silicates. There is so far no detection of planetary 
debris at white dwarfs that has a large C/0 ratio which 
would be indicative of silicon carbide-based minerals. The 
abundances of planetary material found around white dwarfs 
show a large diversity, comparable to, or exceeding that seen 
among different meteorite classes in the solar system. We 
find that the Al/Ca ratio follows a similar trend as observed 
among solar system objects, which suggests that processing 
of proto- and post- planetary material follows similar under- 
lying principles. A particularly interesting pattern is found 
in PG 0843-1-516, where over-abundances of S, Cr, Fe, and Ni 
are suggestive of the accretion of material that underwent 
melting and possibly differentiation. Extending the abun- 
dance studies of metal-polluted white dwarfs both in de- 
tail and number will provide further insight into the diver- 
sity of exo-terrestrial material, and guid e the understand- 
ing of terrestrial exo-p lanet formation l|Bond et al.l I2OI0I : 
ICarter-Bond et aLllioia ). 

^ For completeness, we note that circumstellar high-ionisation 
absorption li nes have also been found around a number of hot 
white dwarfs l lBannister et al1l2003l : Ibickinson et al.ll2012l) . How- 
ever, the origin of the circumstellar material is not clear, and the 
detection of strong C lines suggests a different nature compared 
to the rocky debris found around the stars studied here. 
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